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Abstract. We studied the thermodynamics of the one-dimensional J1-J2 spin-1/2 Heisenberg 
chain for ferromagnetic nearest-neighbor bonds Ji < and frustrating antiferromagnetic next- 
nearest-neighbor bonds J2 > using full diagonalization of finite rings and a second-order 
Green-function formalism. Thereby we focus on J2 < |Ji|/4 where the ground state is still 
ferromagnetic, but the frustration influences the thermodynamic properties. We found that 
their critical indices are not changed by J2. The analysis of the low-temperature behavior of the 
susceptibility x leads to the conclusion that this behavior changes from x^T~^ at J2 < |Ji|/4 
to X T'^^"^ at the quantum-critical point J2 = |>^i|/4. Another effect of the frustration is the 
appearance of an extra low-T maximum in the specific heat Cv{T) for J2 > |Ji|/8, indicating 
its strong influence on the low-energy spectrum. 



Introduction: In low-dimensional frustrated quantum magnets thermal and quantum 
fluctuations strongly influence the low-temperature physics [1,2]. Special attention has been paid 
to one-dimensional (ID) J1-J2 quantum Heisenberg magnets, see Ref. [3] and references therein. 
Recent experimental studies have shown that edge-shared chain cuprates, such as LiVCu04, 
Li(Na)Cu202, Li2ZrCu04, and Li2Cu02 [4~13], represent a family of quantum magnets for 
which the ID J1-J2 Heisenberg model is a good starting point for a theoretical description. The 
above listed compounds are quasi- ID frustrated spin-1/2 magnets with a ferromagnetic (FM) 
nearest-neighbor (NN) in-chain coupling Ji < and an antiferromagnetic (AFM) next-nearest- 
neighbor (NNN) in-chain coupling J2 > 0. 

The model: The Hamiltonian of the ID J1-J2 Heisenberg ferromagnet is given by 



where the first sum runs over the NN bonds and the second sum over the NNN bonds. Henceforth 
we set Ji = —1. For the model ([I]) a quantum critical point at J2 = 0.25 exists where the FM 
ground state (GS) gives way for a singlet GS with spiral correlations for J2 > 0.25 [14-16]. 
For most of the edge-shared chain cuprates J2 is large enough to realize such a spiral GS. 
However, several materials considered as model systems for ID spin-1/2 ferromagnets, such 
as TMCuC[(CH3)4NCuCl3] [17] and p-NPNN (C13H16N3O4) [18], might have also a weak 
frustrating NNN interaction J2 < 0.25. Moreover, recent studies [13] lead to the conclusion 
that Li2Cu02 is a quasi-lD spin-1/2 system with a dominant FM Ji and weak frustrating AFM 




(1) 



J2 ~ 0.2|Ji|. Here we focus on the parameter region J2 < 0.25, i.e. the GS is ferromagnetic. 
Only at J2 = 0.25 the FM GS multiplet is degenerate with a spiral singlet GS [14-16]. On the 
other hand, the frustrating J2 influences the low-energy excitations, in particular, if J2 is close 
to the quantum critical point. Hence, the frustration may have a strong effect on the low-T 
thermodynamics. We mention that previous studies [19, 20] of the thermodynamics of the ID 
J1-J2 model with FM Ji did not consider values of J2 near the quantum critical point J2 ^ 0.25. 

Results: To study the thermodynamic properties we use the full exact diagonalization (ED) 
of finite rings of up to = 22 lattice sites, complemented by data obtained by a spin-rotation- 
invariant second-order Green-function method (RGM) [21-24]. Note that by contrast to ED 
the RGM is limited to values J2 < 0.2 [24] but yields results for 00, that allows the 

calculation of the correlation length by the RGM. Here we will present data for the spin-spin 
correlation functions (SoS.„), the uniform static spin susceptibility x ™d the specific heat C^. 
For the discussion of the correlation length of the model ([1]), see Ref. [24]. For the unfrustrated 
model we will compare our results with available Bethe-ansatz data [25] and transfer-matrix 
renormalization group (TMRG) results [19]. 

The temperature dependence of the spin correlation functions (SoSn) is shown for n = 1 
(NN) and n = 10 for various J2 in Fig. [TJ With increasing frustration the correlation functions 
decrease, where the further-distant correlators decay much stronger than the NN correlator. 
Near and at the quantum critical point the large-distant correlator (SqSio) vanishes already at 
T > 0.05. Interestingly, for J2 = 0.2, 0.24, and 0.25 the correlator (SqSio) changes the sign and 
goes through a minimum. This behavior is not affected by finite-size effects, e.g., the correlators 
(SqSs) for = 16,20 and (SoSg) for = 12,16,20 also change the sign and go through a 
minimum for J2 = 0.2, 0.24, and 0.25. 

Next we discuss the low-temperature properties of the susceptibility x = ^^^h^od{Sz) /dh. 
Due to the FM GS x diverges at T — > 0. Using Bethe-ansatz for J2 = the critical behavior 
has been determined as x oc [25]. Using the RGM, recently it has been confirmed that the 
critical indices for the susceptibility and the correlation length, 7 = 2 and = 1, respectively, 
are not changed by frustration for J2 < 0.25. However, at the quantum critical point J2 = 0.25 
a change of the low-temperature physics is expected [1]. To study that question we consider the 
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Figure 1. Spin correlation function 
(SqSi) (NN) and (SqSio) calculated by 
ED for A^ = 20 sites. 




Figure 2. xT versus vT calculated by ED for 
N = 22 (thick lines - calculated data, thin lines 
- extrapolation to T — > 0, see Eq. ([2]) and text) 
and RGM (circles) as well as Bethe-ansatz data 
(squares) for J2 = from Ref. [25]. The inset 
shows the coefficient uq = lim^^oX^^ versus J2. 
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Figure 3. ED data for the specific heat for 
N = 22. For comparison we show TMRG 
data (squares) from Ref. [19] for J2 = 0. 



Figure 4. Finite-size dependence of the 
specific heat for J2 = 0.22. The inset shows 
CSd and TSo versus J2, see text. 



expansion 

XT^ = yo + yiVT + y2T + 0{T''/^) , (2) 

related to the existence of the FM critical point at T = 0. It has been derived for J2 = 
in Ref. [25]. For the frustrated system ([T]) the coefficients yo, yi, and 2/2 depend on J2. In 
Fig. [2] we plot xT^ versus VT. We find a good agreement of the ED data for with Bethe- 
ansatz results down to quite low temperature T. The RGM results for deviate slightly 
from the Bethe-ansatz results for finite T, but approach the Bethe-ansatz data for T — > 0, 
see also Ref. [22]. The behavior of the leading coefficient yo and the next-order coefficient yi 
can be extracted from the results for by fitting them to Eq. ([2]). For the RGM we use 
data points up to a cut-off temperature T = Tcut = 0.005. To deal with finite-size effects in 
the ED data at very low T, we use the specific heat per site Cy{T), see below, to determine 
that temperature Ted down to which the first four digits of Cy{T) for A'^ = 20 and N = 22 
coincide. Then we fit the ED data in the interval Ted ^ T < Ted + Tcut to Eq. ([2]). Note 
that Ted becomes smaller for increasing J2, we find e.g., Ted = 0.22,0.13,0.09,0.04,0.03,0.02 
at J2 = 0.0,0.1,0.15,0.2,0.24,0.25, respectively. For J2 = we found yo = 1/24 {yo = 0.0418) 
for the RGM (ED), which agrees with the Bethe-ansatz results of Ref. [25]. [Note the different 
definitions of x in oui paper and in Ref. [25].] Including frustration J2 > we observe a linear 
decrease of yo with J2 down to zero at J2 = 0.25 given by 

yo = (1 - 4J2)/24 , (3) 

cf. the inset of Fig. [2j The vanishing of yo at J2 = 0.25 indicates the change of the low-T 
behavior of the physical quantities at the quantum critical point [1]. Indeed, a polynomial fit 
according to yi = ay + byJ2 -|-Cy yields the finite value yi w 0.05 (0.04) for RGM (ED). Hence, 
our data provide evidence for a change of the low-T behavior of x from x T~'^ at J2 < 0.25 to 
X oc T~^/^ at the quantum critical point J2 = 0.25. For a a similar discussion of the correlation 
length ^, see Ref. [24], where it was found that the low-T behavior of ^ changes from ^ oc T~^ 
at J2 < 0.25 to ^ oc T-^/2 at J2 = 0.25. 

In Fig. [3] we present ED results for the specific heat C^. For J2 = we found a broad 
maximum at T ~ 0.332 and a steep decay to zero starting at about T = 0.05 in accord with the 
TMRG [19]. For J2 > 0.125 the specific heat exhibits a minimum located at around T = 0.2, and 
two maxima, namely a high-T maximum at around T = 0.6 and an additional low-T maximum 
at T < 0.1. If J2 approaches J2 = 0.25, a further quite sharp peak at very low T appears. 



that is, however, strongly size dependent, see Fig. [H From Fig. [Hit is obvious that the extra 
low-T finite-size peak behaves monotonously with N. Hence, we have performed a finite-size 
extrapolation to — > oo of the height and the position of the peak in Cy(T) using the 
formula a{N) = Ooo + cli/N"^ + a2/N'^. The extrapolated values and are shown in the 
insets of Fig. [H Obviously, CSd > even near the quantum critical point Ji = 0.25, where 
C£d ~ 0.05. On the other hand, decreases with J2 and becomes very small near J2 = 0.25. 
This behavior suggests that a characteristic steep decay of Cv{T) down to zero starts at very 
low T when approaching J2 = 0.25. 

Summary: We discussed the thermodynamics of frustrated FM spin-1/2 J1-J2 Heisenberg 
chains and found as prominent features (i) a change of the low-T critical behavior at the 
quantum critical point J2 = | Ji|/4, (ii) and an additional low-T maximum in the specific heat 
for |Ji|/4> J2>|Ji|/8. 
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